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ABSTRACT 

The Lya emission has played an important role in detecting high-redshift galaxies, including cur- 
rently the most distant one at redshift z = 8.6. It may also contain important information on the 
origin of these galaxies. Here, we investigate the formation of a typical L* galaxy and its observa- 
tional signatures at the earliest stage, by combining a cosmological hydrodynamic simulation with 
three-dimensional radiative transfer calculations using the newly improved ART 2 code. Our cosmo- 
logical simulation uses the Aquila initial condition which zooms in onto a Milky Way-like halo with 
high resolutions, and our radiative transfer couples multi- wavelength continuum, Lya line, and ioniza- 
tion of hydrogen. We find that the modeled galaxy starts to form at redshift z ~ 24 through efficient 
accretion of cold gas, which produces a strong Lya line with a luminosity of L\, yOL ~ 10 42 erg s _1 as 
early as z ~ 14. The Lya emission appears to trace the cold, dense gas. The lines exhibit asymmetric, 
single-peak profiles, and are shifted to the blue wing, a characteristic feature of gas inflow. Moreover, 
the contribution to the total Lya luminosity by excitation cooling increases with redshift, and it be- 
comes dominant at z > 6. We predict that L* galaxies such as the modeled one may be detected 
at z < 8 by JWST and ALMA with a reasonable integration time. Beyond redshift 12, however, 
only Lya line may be observable by narrow band surveys. Our results suggest that Lya line is one 
of the most powerful tools to detect the first generation of galaxies, and to decipher their formation 
mechanism. 

Subject headings: galaxies: formation - galaxies: evolution - galaxies: high-redshift - radiative transfer 
- line: profiles - hydrodynamics - cosmology: computation 



1. INTRODUCTION 

The quest for the first galaxies formed at the cos- 
mic dawn is a major front ier in both observationa l 
and theoretical cosmology (iBromm fe Yoshidal l2Qllf ). 
Over the past few years, significant progress has 
been made in detecting galaxies at re dshift z > 6, 
using either broad-band colors (e.g., [B ouwens et al I 
ouwens fc Illingwor th 2006; Bouwens et al. 2010, 



201 ll), or narrow - band Lya 
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Malhotra fc Rhoadsl [200l live et all 120061: iStark et al 
2007[lHu et a l. 2010; Leh nert et al . 2010; Vanz ella et al 
2011 fetark et al.ll2011l:IOno et al.ll2012UKashikawa et al 
20111 : iShibuva et al.l 120121 ). In particular, the Lya 
line has played an important role in identifying and 
confirming distant galaxies, the so-called Lya emitters 
(LAEs), including cu rrently the record holder at z = 8.6 
([Lehnert et al.ll2010l ). These remarkable observations in- 
dicate that galaxies formed less than a few hundred mil- 
lion years after the Big Bang. 

Despite the rapidly increasing number of detections, 
the origin and nature of these distant galaxi es, however, 
remain open questions ([Bromm et al .1120091 ). Recently, 
a number of state-of-the-art simulations have started to 



address this issue (e.g.. [Wise fc A bel 2007; Wise et all 
20081 : iWise fc Ab el 2008; IGreif et al.l 1201(1 iWise et al.l 
2012l :IJ eon et al. These studies focused on a halo 

in a small volume (1 Mpc) with high resolutions, and sug- 
gested that the formation of the first galaxies is closely 
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tied to the formation of the first stars (so-called Pop III 
stars) and the feedback from them, and that these galax- 
ies likely consist of second- or third generation of stars 
formed from enriched gas, similar to the present-day 
stars. 

In this work, we explore the physical conditions of early 
galaxy formation on a larger scale. In particular, we fo- 
cus on the gas properties and the Lya emission from it. 
Recent simulations have revealed that a large amount of 
gas penetrate deep inside dark matt e r halos as cold, fil- 

2001 iBirnboim fc Dekell 
Dekel fc Birnboiml l20oT 



amentary streams dKatz et al.l 



2003; Keres et al. 2005, 2009; 
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Ocvirk et al. 2008; Brooks et al. 2009; Dekel et al. 200' 
and iDekel et all (|2009f ) showed that massive galaxies 
at z = 2 — 3 can actively form stars from inflow 
of cold gas. More recently, iDi Matteo et al.l (|201lD 
suggested that massive galaxies at z > 6 can grow 
by cold accretion and evolve with black holes. Such 
streams of cold gas may produce a large number of 
Lya photons via excitation cooling process, and give 
rise to the Lya emission detected in the early galax- 
ies dDiikstra fc L oeb 2009; Fauch er-Gkuere et al.ll20Tol : 
lLatif et a,l-M2Qllb lYaiima et al.ll2012allbl ). 

We combine a multi-scale cosmological hydrodynamic 
simulation with multi- wavelength radiative transfer (RT) 
calculations. The simulation uses the Aquila initial 
condition and follows the formation and evolut i on of 
a Milky Way-size galax y ([Wadepuhl fc Springell 120111 : 
iScannapieco et all |2012[ ) . It covers a large dynamical 
range from a 100 h~ 1 Mpc box down to a ~ 5 h~ 3 Mpc 



zoom-in region, which is ideal to study the gas inflow on 
a large scale. The R T calculations uses th e newly im- 
proved code ART 2 bv lYaiima etHI (|2012a[ ). The ART 2 
code couples multi-wavelength continuum, Lya line, and 
ionization of hydrogen, which is critical to study the Lya 
and multi-band properties of the early galaxies. 

The paper is organized as follows. We describe our cos- 
mological simulation in §2, and the RT calculations in §3. 
In §4, we present the results, which include the Lya prop- 
erties, and detect ability by upcoming missions James 
Webb Space Telescope (JWST) and Atacama Large Mil- 
limeter Array (ALMA). We discuss the implications and 
limitations of our model in §5, and summarize in §6. 

2. GALAXY MODEL 

The cosmological simulation presented here follows the 
formation and evolution of a Milky Way-like galaxy and 
its substructures. The simulation is described in detail 
in (Zhu et al. in prep.), an d it has been employ ed in 
a series of RT calcula tions by lYajima et al.l (j2012af ) and 
lYajima et al.l (j201 2b). Here we briefly highlight some 
important features of the simulation. 

It uses the Aquila i nitial condition, which is t he hy- 
drodynamical version (Wadep uhTfc Springelll201lh of the 
initial condition of the Aquarius Project, the largest 
ever part icle simulation of a M ilky Way-sized dark mat- 
ter halo (Springe l et al.l [20 08). The simulation includes 
dark matter, gasdynamics, star formation, black hole 
growth, and feedback processes, and is performed using 
the parallel, N-body /Smoothed Particle Hydrodynamics 
(SPH) code GAD GET-3, which is an improved ve r sion o f 
that described in lSpringel et al.l (|2001f ): iSpringel (j2005f ). 
GADGET im plements the entropy-conserv ing formula- 
tion of SPH (jSpringel fc Her nauist 2002 ) with adap- 
tive particle smoothing, as in Hernaui st fc Katzl (jl989). 
Radiative cooling and heating processes are calculated 
assuming collisional io nization equilibrium (jKatz et al.l 
1 19961 : iDave et aTlll999[ ). Star formation is modeled in a 
multi-phase ISM , with a rate that follows the S chmidt- 
Kennicutt Law (jSchmidtl 119591 : iKennicuttl Il998h . Feed- 
back from supernovae is captured through a multi-phase 
model of the IS M by an effective equation o f state for 
star- forming gas (Springel & Hernauist 2003). The UV 
background model of Haardt & Madau (1996) is used. 

We also include a model of ga lactic wind driven by 
stellar feedback, as i ntroduced by Springel & Hern quistl 
(|2003D . Similar to iWadepuhl fc Springell (|2011l ). we 
adopt a constant wind velocity of v w ind = 484 km s _1 , a 
mass-loss rate that is twice of the star formation rate, 
and an energy efficiency of unity such that the wind 
carries 100% of the supernova energy. The wind di- 
rection is anisotropical, preferentially perpendicular to 
the galactic disk. This wind model causes an outflow 
of gas, transporting energy, matter and metals out of 
the galactic disk in proportio n to the star formation rate 
(|Springel fc Hernquistl [20031 ) . Our model galaxy shows 
similar SFR, mass of halo, gas and sta rs with those of 
the G3 runs in iScannapieco et~a71 (j2012[ ). 

The whole simulation falls in a periodic box of 
100 /i _1 Mpc on each side with a zoom-in region of a size 
5x5x5 h~ 3 Mpc 3 . The spatial resolution is ~ 250 h~ x 
pc in the zoom-in region. The mass resolution of this 
zoom-in region is 1.8 x 10 6 /i _1 M for dark matter par- 
ticles, 3 x 10 5 /i _1 M for gas, and 1.5 x 10 5 /i _1 M 



for star particles. The cosmological parameters used in 
the simulation are Q m = 0.25, = 0.75, as = 0.9 
and h = 0.73, consistent with t he five-year results of the 
WMAP (jKomatsu et al.ll2009D . The simulation evolves 
from z = 127 to z = 0. 

In the simulation, a galaxy is identified as a group us- 
ing the on-fly friends-of-friends (FOF) group finding al- 
gorithms, which links all baryon particles to their nearest 
dark matter neighbor with a dark matter linking length 
less than 20% of their mean spacing. 

This simulation covers a large dynamical range with 
decent resolutions. It is deal to study the formation of 
the first generation of galaxies and their environments, 
in particular the gas accretion from a large scale. In 
this work, we focus on the most massive galaxy from six 
snapshots at redshift z ~ 6 — 14, representing a Milky 
Way-like galaxy in its infancy. 

3. RADIATIVE TRANSFER 

The RT calculations are performed using the 3D Monte 
Carlo RT code, All-wavelength Radiative Transfer with 
Ad aptive Refinement Tre e (ART 2 ), as recently developed 
by lYajima et al.l (|2012al). ART 2 w as improved over the 
original version of iLi et al.l (j2008[ ) , and features three 
essential modules: continuum emission from X-ray to 
radio, Lya emission from both recombination and col- 
lisional excitation, and ionization of neutral hydrogen. 
The coupling of these three modules, together with an 
adaptive refinement grid, enables a self-consistent and 
accurate calculation of the Lya properties, which depend 
strongly on the UV continuum, ionization structure, and 
dust content of the object. Moreover, it efficiently pro- 
duces mult i- wavelength properties, such as the spectral 
energy distribution and images, for direct comparison 
with multi-band observations. The detail ed implementa- 
tion s of the ART 2 code ar e described in ILi et al.l (|2008f ) 
and lYajima et al.l (j2012af ). Here we focus on the Lya 
calculations and briefly outline the process. 

The Lya emission is generated by two major mecha- 
nisms: recombination of ionizing photons and collisional 
excitation of hydrogen gas. In the recombination process, 
we consider ionization of neutral hydrogen by ionizing 
radiation from stars, active galactic nucleus (AGN), and 
UV background (UVB), as well as by collisions by high- 
temperature gas. The ionized hydrogen atoms then re- 
combine and create Lya photons via the state transition 
2P — >• IS. The Lya emissivity from the recombination is 

4 ec = faOtBhi> a n e nmh (!) 

where an is the case B recombination coefficient, and 
f a is the average number of Lya photons produced 
per case B rec o mbin ation. Here we use ae derived in 
iHui fc Gnedinl (|1997f ). Since the temperature depen- 
dence of f a is not strong, fa = 0. 68 is assumed every- 
where (jOsterbrock fc Ferlandll2006l ) . The product hu a is 
the energy of a Lya photon, 10.2 eV. 

In the process of collisional excitation, high tempera- 
ture electrons can excite the quantum state of hydrogen 
gas by the collision. Due to the large Einstein A coeffi- 
cient, the hydrogen gas can quickly occur de-excitation 
with the Lya emission. The Lya emissivity by the colli- 
sional excitation is estimated by 

e a ° 11 = C Lya n e n Hh (2) 
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where Ch ya is the collisional excitation coefficient, 
fi.vrv = 3.7 x lQ-^expC-fe ^/fcDT- 1 / 2 ergs s" 1 cm 3 
(jOsterbrock fc Ferlandll2006l ). 

Once the ionization structure have been determined, 
we estimate the intrinsic Lya emissivity in each cell by 
the sum of above Lya emissivity, e a = e^ ec + e£° n . 

In RT calculations, dust extinction from the ISM is in- 
cluded. The dust content is estimated according to the 
gas content and metallicity in each cell, which are taken 
from the hydrodynamic simulation. The dust-to-gas ra- 
tio of the MW is used where the metallicity is of Solar 
abundance, and it is linearly interpolated for other metal- 
licity. We u se the stellar population synthesis model of 
GALAXEV (jBruzual fc Charlotll2QQ3l ) to produce intrin- 
sic spectral energy distributions (SEDs) of stars for a grid 
of metallicity and age, and we us e a s imple, br o ken p ower 
law for the AGN ([Li et al.l 120081 ). A lSalpeterl (fl955l ) ini- 
tial mass function is used in our calculations. 

In this work, we apply ART 2 to the selected galaxies 
from the cosmological simulation. In our post-processing 
procedure, we first calculate the RT of ionizing pho- 
tons (A < 912 A) and estimate the ionization fraction 
of the ISM. The resulting ionization structure is then 
used to run the Lya RT to derive the emissivity, fol- 
lowed by the calculation of non-ionizing continuum pho- 
tons (A > 912 A) in each cell. Our fiducial run is done 
with Afph = 10 5 photon packets for each ionizing, Lya, 
and non-ionizing compone nts, which was demo nstrated 
to show good convergence ([Yajima et al.l 12012 al ibi). The 
highest refinement of the adaptive grid corresponds to a 
cell size comparable to the spatial resolution of 250 pc in 
comoving coordinate of the hydrodynamic simulation. 

4. RESULTS 

The formation history of a Milky Way-like galaxy, 
its multi-band properties, and its connection to high- 
redshift Lya emitters, have b een presented in a s eries 
of papers (Zhu et al. in prep., lYaiima et aT1l2M2l bllcT>. 
Here we focus on the galaxy formation at its earliest evo- 
lutionary stage and the Lya properties from z ~ 6 — 14. 

4.1. The Accretion of Cold Gas 

The modeled MW galaxy starts to form at z ~ 24. 
Figure [T] shows the distribution of gas density, gas tem- 
perature, and stellar density of the MW main progenitor 
from redshift z ~ 14 to z ~ 6. The gas follows the dis- 
tribution of dark matter and exhibits filamentary struc- 
tures. At z > 6, the gas is predominantly cold, with a 
mean temperature of ~ 10 4 K. Stars form out of such 
cold gas, so they also distribute along the filaments. 

The star formation at z > 6 is fueled by efficient accre- 
tion of cold g as, as demonstrated in Figure [2] The gas 
accretion rate is defined as the inflow rate of gas within 
the virial radius of the modeled galaxy. It peaks around 
10 4 K in all cases. At a later time, feedback from both 
stars and accreting BHs heats up the gas. Also, the gas 
can be heated by gravitational shocks during the infall. 
Therefore the accretion includes hot gas as well. The in- 
flow gas falls along the filaments toward the intersection, 
the highest density peak where the first galaxy in the 
simulated volume forms. 

Figure [3] shows the star formation history of the MW. 
The star formation rate (SFR) increases steadily from 



- 3 x 10 -3 M yr _1 at z ~ 24 to ~ 15 M yr _1 at 
z ~ 8.5, and it peaks at ~ 62 M yr _1 at z = 5.2, owing 
to a merger of gas-rich protogalaxies. The galaxy mass 
increases rapidly during this cold accretion phase. By 
z ~ 8.5, it reaches a total mass of ~ 5.6 x 10 10 M , and 
a stellar mass of ~ 6 x 10 9 M . 

4.2. The Lya Properties 

The Lya emission traces the gas distribution, as 
shown in Figure |U The surface brightness rises above 
10 -20 ergs s _1 cm -2 arcsec -2 at z < 14. At high red- 
shift z > 10, the galaxy is small, and the Lya emission 
is faint and confined to the central high-density region. 
As the galaxy grows in mass and size, the Lya emission 
becomes stronger and more extended and irregular due 
to mergers and gas infall along the filaments of the main 
halo. 

Figure [5] shows the Lya properties of the MW galaxy 
at from z ~ 14 to z ~ 6, including the emergent Lya 
luminosity (^Lya), equivalent width of Lya line in rest 
frame, and photon escape fraction of Lya and UV con- 
tinuum (1300 - 1600 A). For comparison with the star 
formation activity, the SFR of the galaxy at correspond- 
ing redshift is also shown. 

During this early growth phase, the SFR of the galaxy 
increases from ~ 1 M yr _1 at z ~ 14.0 to ~ 31 M yr _1 
at z ~ 6, owing to abundant supply of cold gas from 
infall and merging of gas-rich mini halos. The result- 
ing emergent Lya luminosity shows a similar trend, in- 
creasing from ~ 1.6 x 10 42 ergs _1 at z ~ 14.0 to ~ 
5.5 x 10 42 erg s _1 at z ~ 6. If we consider only the recom- 
bination process with the assumption of Li jya /Ln a = 8.7 
(in which the H a is a tracer of star formation), the the 
intrinsic Lya luminosity should be linearly proportional 
to SFR, Lt yfy (er g s" 1 ) = 1.1 x 10 42 x SFR (M yr" 1 ) 
(jKennicu tt 1998). However, the evolution of L\, ya in 
Figure [5] differs from the SFR history. This is due to the 
contribution from excitation cooling to the Lya emission 
as we will discuss later, and dust absorption of the Lya 
photons. In particular, at z > 10, the L\, ya increases 
with redshift, in opposite direction from the SFR, as a 
result of high collisional excitation and high / esc . 

The lower-left panel of Figure [5] shows the photon 
escape fraction of Lya, and the UV continuum 

f^, where is calculated at A res t = 1300 - 1600 A. 
The fe/ c a of the modeled galaxy falls in the range of 
0.49 — 0.81, and it increases with redshift. In ou r model, 
the d ust is produced by type-II supernovae (|Li et al.l 
120081 ) . The dust amount increases as star formation rises 
from z ~ 14 to z ~ 6, and hence it efficiently absorbs the 
Lya and UV continuum photons, resulting decreasing es- 
cape fraction. However even at z = 8.5, about 40% of 
Lya photons are absorbed by dust. This is due to the 
fact that, in the early phase, galaxies are gas rich and 
compact, the gas and dust are highly concentrated in 
the galaxies, resulting in effective absorption of the the 
Lya and UV photons by the dust. 

The resulting Lya equivalent width (EW) is shown in 
the lower-right panel of Figure [5] The EW is defined 
by the ratio between the Lya flux and the UV flux den- 
sity fxjy in rest frame, where the mean flux density of 
A = 1300 — 1600 A is used. The modeled galaxy has 
EW > 20 A at these redshifts, and is therefore classi- 
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Figure 1. The distribution of gas density (left column), gas temperature (middle column), and stellar density (right column), of the MW 
galaxy at z ~ 14, 10.4,8.5, and 6.2, respectively. The box size is 1 Mpc in comoving scale. The temperature of the gas is in Kelvin in log 
scale, as indicated in the color bar. 
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Figure 2. The gas accretion rate by the MW galaxy as a function 
of gas temperature at different redshift. 
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Figure 3. The growth history of the MW galaxy illustrated by 
the star formation rate (top panel), and by the accumulated mass 
(bottom panel), in which filled circle represents the total mass, 
while the open circle represents the stellar mass. 



fied as a Lya emitter (LAE) (e.g.. IGronwall et alJ[2007l ). 
The EW increases with redshift, from ~ 93 A at z ~ 6 
to rsj 2300 A at z ~ 14.0. Such trend is in broad agree- 
ment with observations that galaxies at higher redshifts 
appear to have highe r EWs than their co u nterparts at 
lower redshifts (e.g., iGronwall et al.l [20071 : lOuchi et al.l 
120081 ). This is because the contribution from excitation 
Lya cooling becomes large with increasing redshift, as 



shown in Figure [6] in the next section, which boosts the 
EW significantly (|Yajima et alJl2Q12aH bh . 

The currently most distant galaxy detected from Lya 
emission line at z = 8.6, UDFy- 38135539, has a Lya 
luminosity of - 5.5 x 10 42 ergss" 1 (jLehnert et al.ll201Q[ ). 
It is somewhat brighter than our modeled galaxy at the 
same redshift, indicating that the host galaxy may be 
more massive than the one in our simulation. 

4.3. Contribution of Excitation Cooling to Lya 
Emission 

As mentioned in § O the Lya emission is generally pro- 
duced by the recombination of ionizing photons and the 
collisional excitation of hydrogen gas. In our cosmologi- 
cal simulation, galaxy evolution is accompanied by cold, 
filamentary gas streams with temperature T ~ 10 4-5 K, 
which penetra te deep inside the da rk matter halos (Zhu 
et al. in prep, lYajima et al.l 12012 rih a phenomenon also 



reported by other groups (Katz et al. 2003; K eres et al 
2005, 2009; Birnboim & Dekel 2003; Dekel & Birnboim 
2006; lOcvirk et al"1l2008t iBrooks et aD2009l : iDekel et al 
20091 ). Such cold gas can effic iently produce the ex- 
citation Lya cooling photons (pD iikstra & Loebl 120091 : 
iFaucher-Gkuere et al.l 120091 : IGoerdt et al.l 12010( 1. At 
higher redshift, galaxies experience more merging events 
and accrete more cold gas efficiently, which results 
in stronger Lya emission from excita tion cooling, and 
higher Lya EWs (jYaiima et al.l 120 12al lbh. 

As shown in Figure [6j the fraction of excitation cooling 
Lya to the total intrinsic Lya luminosity increases from 
~ 65% at z ~ 6 to ~ 88% at z ~ 14. Such extremely 
high excitation Lya cooling produces the extremely high 
Lya EWs seen in Figure [51 

4.4. The Lya Line Profile 

The resulting Lya line profiles of the modeled MW 
galaxy from z ~ 14 to z ~ 6 are shown in Figure 
The frequency of the intrinsic Lya photon is sampled 
from a Maxwellian distribution with the gas tempera- 
ture at the emission location in the rest flame of the 
gas. All sources show asymmetric profiles with a sin- 
gle peak or weak double peaks. More interestingly, 
most profiles are shifted to the shorter (bluer) wave- 
lengths. This is a characterist ic feature of gas inflow 
(Zheng & Mirald aTEscudel l2002f ). Indeed, as shown in 
Figure [HI a significant fraction of the gas shows a large 



infalling velocity V r 



-100 to 



-200 km s , even 



though our simula tion includes feedbac k of stellar wind 
similar to that of iSpringel et al.l (|2005[ ). In particular, 
the gas in the galaxy from z ~ 14 — 10 is dominated by 
inflow motion, which explains the significant blue shift of 
the profiles in Figure [7] (top panel). At redshift z < 8.5, 
the gas exhibits outflow as well, and has a larger velocity 
distribution —250 < V r < 200 km s _1 , which results in 
an extended profile to both blue and red wings. 

We note that the Lya line profile may be suppressed 
and changed by the intergalactic medium (IGM) (e.g. , 
ISantosI l200l iDiikstra et al.l 120071 : iZheng et al.l 120101 : 
lLaursen et al.ll2011h . because the IGM effectively scat- 
ters the Lya photons at the line center and at shorter 
wave lengths by the Hubble flow (e.g., lLaursen et al.l 
120111 ). As a result, the inflow feature in our profiles may 
disappear and the shape may become an asymmetric sin- 
gle peak with only photons at red wing. However, while 
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Figure 4. The evolution of the Lyo; surface brightness of the MW galaxy with redshift, at z ~ 14, 10.4, 8.5, and 6.2, respectively. The 
box size is 1 Mpc in comoving scale. The color indicates the Lya surface brightness in log scale in units of ergs s _1 cm -2 arcsec -2 . 



asymmetric line profiles with an extended red wing are 
commonly seen in high-redshift LAEs, there appears to 
be some profiles in the z > 6 observations that have 
complex features including double peaks a nd extended 
blue wing, similar to what we see here (e.g.jOuchi et al.l 
[20101 : IHu et all 1201(1 iKashikawa etalll2Qilh , The ob- 
served line of UDFy-38135539, the currently highest red- 
shift galaxy detected at z = 8.6, is not resolv ed and thus 
has a Gaussian profile (Leh nert et al.l [2Q10h . More ob- 
servations of high-resolution Lya line profiles of high- 
redshift LAEs are needed to test our model and verify 
our predictions. 

4.5. Detectability of Progenitors of Local L* Galaxies 

The corresponding multi-wavelength SED of the MW 
galaxy at different redshift are shown in Figure [9j The 
shape of the SED changes significantly from z = 14.0 to 
z = 6.2, as a result of the changes in radiation from stars, 
absorption of ionizing photons by gas and dust, and re- 
emission by dust in the infrared. The galaxy shows Lya 
line in all cases, and the deep decline of UV continuum 



at z > 8.5 is caused by strong absorption of ionizing 
photons from the dense gas. 

A major science goal of the two forthcoming telescopes, 
ALMA and JWST, is to detect the first galaxies. In order 
to predict the detectability of the infancy of a local L* 
galaxy, we contrast the SEDs with some detection limits 
of these two facilities in Figure [9] Our calculations show 
that the flux at 850 /am in the observed frame of the 
model galaxy ranges from ~ 7.9 x 10 -5 mJy at z = 14.0 
to ~ 4.7 x 10 -2 mJy at z = 6.2. With an array of 50 
antennas and an integration of 10 hours, ALMA may 
be able to detect such galaxies at z < 8.5 with a 3 <r 
significance. However, since galaxies do not have a lot 
of young stars and much dust at z > 10, observations 
in continuum by ALMA becomes more difficult, and it 
would need tens of hours of integration time. In contrast, 
JWST appears to be more powerful to detect the earliest 
galaxies as the one we model here, because it can detect 
the UV continuum in rest frame up to z ~ 10. The 
Lya emission is strong even at z ~ 12, which may be 
observable by narrow band surveys with similar detection 
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Figure 5. The Lya properties of the modeled galaxy from z ~ 14 
to 2 ~ 6, including, in clockwise direction, star formation rate, 
emergent Lya luminosity, equivalent width of Lya line in rest 
frame, and photon escape fraction of Lya (filled circles) and UV 
continuum (1300 - 1600 A, open circles). 
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Figure 7. The Lya line profile of the modeled galaxy at different 
redshifts. The dotted- and solid lines are the intrinsic and emergent 
Lya profile, respectively. 
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Figure 6. The fraction of excitation cooling Lya to the total 
intrinsic Lya luminosity as a function of redshift from z ~ 14 to 
z ~ 6. 
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Figure 8. The probability distribution function of the gas mass of 
the neutral hydrogen in the galaxy as a function of radial velocity. 
The velocity is estimated from the center of mass of the galaxy in 
the radial direction. 



limit as the JWST. For example, if we use the F164N 
filter of JWST, the Lya flux of our galaxy at z ~ 12 will 
be detected at 3 a with an integration of ~ 20 hours. 

We note that in the above estimation, IGM absorp- 
tion and transmission were not taken into account. 
The IGM can significantly suppress the Lya flux, and 



the transmission hi ghly depends on viewing angle (e.g., 
lLaursen et al.l 120111) by inhomogeneous ionization struc- 
ture in IGM (e.g.. I Abel et al.l 120071: lYoshida et1Hl200ft 
Ueeson-Daniel et al.l I2012L which make the detection 
more difficult. Of course, the galaxies we present here 
represent progenitor of a local L* galaxy such as the 
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Figure 9. The spectral energy distribution of the MW galaxy 
at different redshift and its detect ability with JWST and ALMA. 
The open squares indicate the 10 a detection limits of JWST and 
ALMA with 16 antennas at an integration time of 1 and 10 hours 
(from top to bottom), while the open triangles indicate the 3 a 
detection limits of JWST and ALMA with 50 antennas. 



Milky Way. Galaxies formed in highly overde nse regions 
are likely much more massive ([Li et alj|2007[ h and may 
be more easily detected by both ALMA and JWST (Li 
et al, in preparation). 

5. DISCUSSIONS 

The gas inflow feature is present in our simulation with 
outflow from stellar feedback. In order to probe the effect 
of stellar wind on the gas inflow, we also performed the 
simulation with a pure thermal feedback model, in which 
the feedback from supernovae is only in thermal energy. 
In such a model, some fraction of thermal energy can 
quickly escape as cooling radiation before conversion to 
kinetic energy. As a result, gas outflow does not occur 
efficiently. 

The resulting Lya line profiles and the probability dis- 
tribution function of the neutral gas mass are shown in 
Figure [10] and Figure [TTJ respectively. Without a strong 
outflow, the Lya line profiles show more pronounced blue 
wing. 

The galaxy in our simulation resides in a low over- 
density region and it represents those that would evolve 
into present-day L* galaxies such as the MW, so it may 
not be the very first one formed in the universe. It is 
believed that the most massive halos in the highly over- 
dense ij^ejyxDms^^ stars form 
(e.g lAbel et alJl2 QQ2: Bro mm fc LarsoiJ2QQ4 : iGao et al.l 
120071 : lYoshida et aTll2008[ ). These may also be the for- 
mation sites of the very first galaxies, owing to feedback 
and chemical enrichm ent from the PopIII stars, as well a s 
abundant gas supply (|Li et al.l 2007: B romm et al.ll2009l ). 

One of the major limitations of our model is that the 
cosmological simulation does not have sufficient resolu- 
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Figure 

include 



10. Same as in Figured but here the simulation does not 
the wind model from stellar feedback. 
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Figure 11. Same as Figure [8] but here the simulation does not 
include the wind model from stellar feedback. The probability 
distribution function of the gas mass of the neutral hydrogen in the 
galaxy as a function of radial velocity. The velocity is estimated 
from the center of mass of the galaxy in the radial direction. 



tions to follow the formation and evolution of individual 
stars. Instead, star formation is modeled using a "sub- 
grid" recipe based on the observed Schmidt-Kennicutt 
Law (Kennicutt 1998). Gas particles are converted into 
stars once it is coole d below 10 4 K and the d ensity 
is above a threshold (Springel & Hernquist 2003). Al- 
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though this treatment is rather simplistic, it neverthe- 
less gives a global star formation history close to what is 
believed of the MW galaxy. 

Another major limitation is that our current RT calcu- 
lations do not include the propagation and scattering of 
Lya and ionizing photons in the IGM. We make the pre- 
diction that galaxies with inflow of cold gas would result 
in asymmetric, blue-shifted Lya line profiles. However, 
as discussed earlier, the absorption by IGM may change 
the the profile to one with extended red wing. We will 
study this issue in more detail in future work that in- 
cludes the radiative transfer of Lya and ionizing photons 
in the IGM. 

6. SUMMARY 

In this work, we have investigated the formation of 
a typical, nearby L* galaxy such as the MW, and its 
Lya properties at the earliest evolutionary stage. We 
combine a cosmological hydrodynamic simulation, which 
uses the Aquila initial condition and focuses on a MW- 
like galaxy, with three-dimensional radiative transfer cal- 
culations using the improved ART 2 code, which couples 
multi- wavelength continuum, Lya line, and ionization of 
hydrogen. 

We find that the modeled MW galaxy forms from ef- 
ficient accretion of cold gas early on, which sustains a 
high star formation rate from z ~ 14 — 6. The cold ac- 
cretion produces strong Lya emission via collisional exci- 
tation, which has a luminosity from ^ 1.6 x 10 42 erg s _1 
at z ~ 14 to - 5.5 x 10 42 erg s" 1 at z ~ 6. The es- 
cape fraction of Lya photons increases from ~ 0.49 at 
z ~ 6 to ~ 0.81 at z ~ 14, due to less dust content at 
higher redshift. The EWs of the Lya lines increases with 
redshift, from - 93 A at z ~ 6 to - 2300 A at z ~ 14. 
Such high EWs may be due to significant contribution to 
Lya emission by excitation cooling, which dominates at 
high redshift. The resulting Lya lines exhibit asymmet- 
ric, mostly single-peak profiles shifted to the blue wing, 
a characteristic feature of inflow. 

Furthermore, we demonstrate that progenitors of local 
L* galaxies such as the modeled one may be detected at 
z < 8 by JWST and ALMA with a reasonable integration 
time. At higher redshift z > 12, however, only Lya line 
may be observable by narrow band surveys with similar 
detection limit as JWST. 

Our results suggest that Lya line may be used to probe 
the formation and evolution, and gas properties of dis- 
tant galaxies. It is perhaps one of the most powerful tools 
to detect the first generation of galaxies in the coming 
decade. 
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